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Abstract 

We derive the triply differential spectrum for the inclusive rare decay B — ► X s l + £~ in the shape 
function region, in which X s is jet-like with m\ < m;,AQCD- Experimental cuts make this a relevant 
region. The perturbative and non-perturbative parts of the matrix elements can be defined with 
the Soft-Collinear Effective Theory, which is used to incorporate a s corrections consistently. We 
prove that, with a suitable power counting for the dilepton invariant mass, the same universal jet 
and shape functions appear as in B — > X s ^f and B — > X u tv decays. Parts of the usual a s (mb) 
corrections go into the jet function at a lower scale, and parts go into the non-perturbative shape 
function. For B — ► X s £ + £~ , the perturbative series in a s are of a different character above and 
below [i = m&. We introduce a "split matching" method that allows the series in these regions to 
be treated independently. 
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I. INTRODUCTION 



The B meson is particularly suitable for probing QCD and flavor physics in the Standard 
Model, since the large mass of the b quark relative to Aq CD provides a useful expansion pa- 
rameter, ~ 0.1. The study of inclusive B decays circumvents the need for precision 
hadronic form factors, while still allowing model-independent predictions. Rare inclusive de- 
cays, which involve flavor- changing neutral currents (FCNCs), not only allow measurements 
of CKM matrix elements, in particular V ts and V t d, but are also highly sensitive to new 
physics, since they do not occur at tree level in the Standard Model. 

Among the inclusive rare B decays, the radiative process B — > X s 7 has received the most 
attention, having been measured first by CLEO (T| and subsequently by other experiments 
0, |3|, [|, |5| . These measurements have provided significant constraints on extensions to the 
Standard Model. The decay B — > X s £ + £~ is complementary to, and more complicated than, 
B — > Xg'y. Its potential for revealing information beyond that supplied by the radiative 
decay is due to the presence of two extra operators in the effective electroweak Hamiltonian 
and the availability of additional kinematical variables, such as the dilepton invariant-mass 
spectrum and the forward-backward asymmetry. Belle and BABAR have already made initial 
measurements of this dilepton process [§, [/J. 

Provided that one makes suitable phase-space cuts to avoid cc resonances, B — > X s £ + £~ 
is dominated by the quark-level process, which was calculated in Ref. . Owing to the dis- 
parate scales, m b <C mw, one encounters large logarithms of the form a™ (m^) \og n (mb/mw) 
(leading log [LL]), a™ +1 (m b ) log 1 (mi, / m w ) (next-to-leading log [NLL]), which should 



be summed. The NLL calculations were completed in Refs. [[10], [TT|| , and the NNLL analysis, 
although technically not fully complete, is at a level that the scale uncertainties have been 



substantially reduced, after the combined efforts of a number of groups [12, 13, [L4L [L5L 16 



Non-perturbative corrections to the quark-level result can also be calculated by means 
of a local operator product expansion (OPE) fl?|| , with non-perturbative matrix elements 
defined with the help of the Heavy Quark Effective Theory (HQET) |l8fl . As is the case 
for B — > Xg'y and B — > X u £u, there are no OiXjm^) corrections. The 0(l/m 2 ) corrections 
and OPE were considered in Ref. []TU[ and subsequently corrected in Ref. PD| . The 0(1/ mp 



corrections were computed in Ref. [21]. There are also non-perturbative contributions arising 
from the cc intermediate states. The largest cc resonances, i.e. the J/ip and if)', can be 
removed by suitable cuts in the dilepton mass spectrum. It is generally believed that the 
operator product expansion holds for the computation of the dilepton invariant mass as 
long as one avoids the region with the first two narrow resonances, although no complete 
proof of this (for the full operator basis) has been given. A picture for the structure of 
resonances can be obtained using the model of Kriiger and Sehgal ||22|| , which estimates 
factorizable contributions based on a dispersion relation and experimental data on cr(e + e~ — > 
cc + hadrons). Non-factorizable effects have been estimated in a model-independent way by 



means of an expansion in l/m c [23], which is valid only away from the resonances. 

Staying away from the resonance regions in the dilepton mass spectrum leaves two 
perturbative windows, the low- and high-g 2 regions, corresponding to q 2 < 6 GeV 2 and 



■ 2 > 14.4 GeV 2 respectively. These have complementary advantages and disadvantages |16 



For example, the latter has significant 1/rrib corrections but negligible scale and charm-mass 
dependence, whereas the former has small 1/rrib corrections but non- negligible scale and 
charm-mass dependence. The low-g 2 region has a high rate compared to the high-g 2 region 
and so experimental spectra will become precise for this region first. However, at low q 2 
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an additional cut is required, making measurements less inclusive. In particular, a hadronic 
invariant-mass cut is imposed in order to eliminate the combinatorial background, which 
includes the semileptonic decay b — > c(— > se + v)e^v = b — > se + e~ + missing energy. The 
latest analyses from BABAR and Belle impose cuts of m x < 1.8 GeV and < 2.0 GeV 
respectively || 0, ||, which in the 5-meson rest frame correspond to g° > 2.3 GeV and put 
the decay rate in the so-called shape function region pif . This cut dependence has so far 
been analyzed only in the Fermi-motion model [EH . 



Existing calculations for B — ► X s £ + £~ are based on a local operator product expan- 
sion in Aqcd/^6- When m 2 x < m^A ~ (2 GeV) 2 , this operator product expansion breaks 
down, and, instead of depending on non-perturbative parameters (Ai, A2, • • •) that are ma- 
trix elements of local operators, the decay rates depend on non-perturbative functions. 
Furthermore, in this region the standard perturbative a s corrections to the partonic process 
b — > s£ + £~ do not apply, since some of these corrections become non-perturbative. Thus, 
even at leading order there does not exist in the literature a model-independent computation 
of the B — ► X s £ + £~ decay rate that can be compared directly with the data at low q 2 . 

Here we study B — > X s £ + £~~ (£ = e, p) in the shape function region for the first time. 
The relevant scales are m 2 ^ 3> m\ 3> mbAQCD 3> Aq CD . In this paper we derive the proper 
theoretical expression for the leading-order triply differential decay rate, which incorporates 
non-perturbative effects that appear at this order and a correct treatment of the perturbative 
corrections at each of the scales. Using the Soft-Collinear Effective Theory (SCET) |26|, 
|2"T| , p8| , |29|| we prove that the non-perturbative dynamics governing the measurable low-g 2 
spectra in B — > X s £ + £~ is determined by the same universal shape function as in endpoint 
B — > X u £v and B — > Xpj decays. We also prove that the decay rate can be split into a 
product of scale-invariant terms, capturing physics at scales above and below m^. We show 
that this procedure, which we call "split matching", can be used to deal with a tension 
between the perturbative corrections that come from these two regions. Implications for 
relating the B — > X s £ + £~ measurements with the mx cut to the Wilson coefficients are 
presented in a companion publication |5D] . 

In the shape function region, the set of outgoing hadronic states becomes jet-like and the 
relevant degrees of freedom are collinear and ultrasoft modes. This is why the appropriate 
theoretical method is SCET. The endpoint region has been the focus of much work in the 



context of B -> X sl and B -> XJv (see e.g. Refs. f2|, ||, [H], ||, g|, ||, H, H, |37|, g|, [3S 



40, 41, 42, 43]). In B — > X u £v this is because of the cuts used to eliminate the dominant 
b — ► c background. In B — > X s j, it is known that cuts with q° •> 2.1 GeV put us in the 
shape function region. 1 

In the small-g 2 region of B — ► X s £ + £~~ with q° > 2.3 GeV, shape-function effects also 
dominate rather than the expansion in local operators. To see this, we note that the mx cut 
causes 2msEx = m 2 B +m 2 x — q 2 ^> m\. Decomposing 2Ex = p x +Px w hh m\ = p x Pxi we 
see that the X s is jet-like with p x S> p x , and the restricted sum over states in the X s causes 
the non-perturbative shape functions to become important. For the experimental cuts on 
q 2 and m x , values for p x are shown in Fig. [I]. It should be clear from this figure that the 



1 In Rcf. Q it was pointed out that even a cut of E n > Eq = 1.8 GeV, corresponding to mx ^ 3 GeV, might 
not guarantee that a theoretical description in terms of the local OPE is sufficient, owing to sensitivity 
to the scale A = mb ~ "2Eq in power and perturbative corrections. Using a multi-scale OPE with an 
expansion in A/ A allows the shape function and local OPE regions to be connected [[33], [54], p| . 
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FIG. 1: The kinematic range for p x and p x given the experimental cuts of q 2 < 6GeV 2 and 
m x <2.0GeVfor B^X s £+£~. 

measurable spectrum is dominated by decays for which p x 3> p x . 

To compute B — > X s £ + £~ in the shape function region with renormalizat ion-group evo- 
lution requires the following steps: 

i) matching the Standard Model at /i ~ m\y on to H\y, 

ii) running H w to /i ~ m^, 

iii) matching at fi ~ on to operators in SCET, 

iv) running in SCET to /i ~ ^/m^Ji, 

v) computation of the imaginary part of forward-scattering time-ordered products in 
SCET at fi ~ y/m^A. This leads to a separation of scales in a factorization theorem, 
which at LO takes the form 2 

d 3 T^=H JdkjW(k)f(°\k), 

with perturbative H and and the LO non-perturbative shape function /(°), 

vi) evolution of the shape function /(°) from Aqcd up to fi ~ a/"^Aqcd • 

For the shape-function decay rate, steps i-ii) are the same as the local OPE results for B — > 
X s £ + £~ . Furthermore, based on the structure of leading-order SCET operators that we find 
for B —>■ X s £ + £~ , we demonstrate that results for other inclusive endpoint analyses can be 



2 Note that the operator product expansion used here occurs at [i ~ rather than at to 2 , as in the 

standard local OPE. 
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used in steps iv) and vi) p6|, [27], 0. 3 Because of this our computations focussed on steps iii) 



and v). In step iii) we show how to implement the split-matching procedure to formulate the 
perturbative corrections, which we elaborate on below. In step v) we derive a factorization 
theorem for B — > X s £ + £~ . This includes computing the hard coefficient functions H at NLL 
order and formulating the structure of these terms to all orders in a s . It also includes a 
derivation of formulae for the decay rate and forward-backward asymmetry that properly 
take into account the effect of the current experimental cuts and the perturbative and non- 
perturbative corrections. 

At leading order in the power expansion the result of steps i)-vi) takes the schematic form 

d 3 T (0) = £{^w)Uw{^w^o)B{^)U H {^Q^i)J{^i)Us{^i^k)f {0 \^K), 

fi w ~ m w , fi ~ m b , fii ~ (m fe A) 1/2 , /x A ~ 1 GeV , (1) 

where £ , B and J represent matching at various scales, and Uw, Uh and Us represent the 
running between these scales. Eq. (JJ) shows only the scale dependence explicitly, not the 
kinematic dependences or the convolutions between J ', Us, and f(°\ which we describe later 
on. 

In a standard application of renormalization-group-improved perturbation theory (LL, 
NLL, NNLL, etc.), the results at each stage of matching and running are tied together, as 
depicted in Eq. ([]]). Usually this would not be a problem, but for B — > X s £ + £~ the nature 
of the perturbative expansion above and below fi ~ rrib is different. Above fi ~ mj the 
series of (a s ln) fc terms are of the traditional form, with a basis of ~ 10 operators (including 
four-quark operators), whose mixing is crucial. Below fi ~ rrib we demonstrate that the 
evolution is universal (to all orders in a s ) for the leading-order operators, but there are 
Sudakov double logarithms of the ratios of scales, which give a more complicated series. It 
turns out to be convenient to decouple these two stages of resummation so that one can 
consider working to different orders in the a s expansion above and below fi = rrib- There is a 
simple reason why this decoupling is important: for ii > m& the power counting and running 
are for currents in the electroweak Hamiltonian and dictate treating Cg ~ l/at s with C-j ~ 1 
and Cio ~ 1. However, at /i = m;, the coefficients Cg and C\g are numerically comparable. 
For fi < nib in the shape function region we must organize the power counting and running 
for time-ordered products of currents in SCET rather than amplitudes, and it would be 
vexing to have to include terms oc C g 2 to 0(a 2 s ) before including the C 2 and C 7 2 terms at 
order 0(oP s ). Thus, once we are below the scale m b , a counting with Cg ~ C w ~ C 7 ~ 1 is 
more appropriate. 

To decouple these two regions for B — > X s £ + £~ decays we make use of two facts: i) for 
fi > rrib the operator Oio involves a conserved current and has no operators mixing into it, 
so it does not have an anomalous dimension, and ii) for ii < rrib all LO biquark operators 



in the Soft-Collinear Effective Theory have the same anomalous dimension ||27|| . We shall 
show that the operators for B — > X s l + l~ are related to these biquark operators. These 
properties ensure that we can separate the perturbative treatments in these two regions at 



In step iv) we can run the hard functions down using results from Refs. [|6H3- In step vi) we can run 
the shape function up to the intermediate scale using the simple result from Ref. |34| ]. An equally valid 
option would be to evolve the perturbative parts of the rate down to a scale /i ~ 1 GeV, as considered 
earlier El IH, Ei El . 
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any order in perturbation theory. This is done by introducing two matching scales, [1q ~ mj 
and fib — trib- The two aforementioned facts allow us to write 

Uw(nw, ^o)<6(/x )£M/i ,/ii) = U w (fi w , fi )B(iJ,o, Vb)U H (lJ>b, Vi) 

= U w (fi W , fio)B 1 (fi )B 2 {fi b )U H (fi b , fii) , (2) 

with well-defined B 1 and B 2 . We define B 2 (fi b ) by using the matching for the operator O w 
and extend this to find B 2 matching coefficients for the other operators using property ii) 
above. The remaining contributions match on to B\. Diagrams which are related to the 
anomalous dimension for fi > uib end up being matched at the scale fio on to Bi, while those 
related to anomalous dimensions for fj, < irtb are matched at a different scale, fib, on to B 2 . 
This leaves 



d 3 r (0) 



£{fi w )U w {n w , fi )Bi{fi ) B 2 (fib)U H (fib, Vi)J(fii)U s (ni, /i A ) / (0) (/x A ) 



(3) 



which is the product of two pieces that are separately //-independent. We refer to this 
procedure as "split matching" because formally we match diagrams at two scales rather 
than at a single scale. The two matching /i's are "split" because they are parametrically 
similar in the power-counting sense. 

We organize the remainder of our paper as follows. We begin by using split matching to 
determine the hard matching functions, B = BiB 2 , for B — > X s £ + £~ in SCET; this is one of 
the main points of our paper. It is discussed in Sec. [11 A| at leading power and one-loop order 
(including both bottom-, charm-, and light-quark loops and other virtual corrections). The 
extension to higher orders is also illustrated. Steps i) and ii) are summarized in Sec. [II A| , 
together with Appendix In Sec. [11 B| we discuss the running for step iv) and give a brief 
derivation of why the anomalous dimension is independent of the Dirac structure to all 
orders in a s . In Sec. [11 C , we discuss the basic ingredients for the triply differential decay 
rate and the forward-backward asymmetry in terms of hadronic tensors. A second main 
point of our paper is the SCET matrix-element computation for B — > X s £ + £~, step v), 



which is performed in Sec. |il U| . In Sec. [11 E| we review the running for the shape function, 
step vi). In Sec. [TlJ we present our final results for the differential decay rates at leading 
order in the power expansion, including all the ingredients from Sec. [TT] and incorporating 
the relevant experimental cuts. The triply differential spectrum and doubly differential 
spectra are derived in subsections [IljA-D. Readers interested only in our final results may 



skip directly to section |T|. We compare numerical results for matching coefficients at m?, 
with terms in the local OPE in Sec. |IIE[ In Appendix || we briefly comment on how our 
analysis will change if we assume a parametrically small dilepton invariant mass, q 2 ~ A 2 , 
rather than the scaling q 2 ~ A used in the body of the paper. (For the case q 2 ~ A 2 , the 
rate for B — > X s £ + £~ would not be determined by a factorization theorem with the same 
structure as for B — > X u £p.) 



II. ANALYSIS IN THE SHAPE FUNCTION REGION 
A. Matching on to SCET 

We begin by reviewing the form of the electroweak Hamiltonian obtained after evolution 
down to the scale /i ~ rrib, and then perform the leading-order matching of this Hamiltonian 
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on to operators in SCET. For the treatment of 75 we use the NDR scheme throughout. 
Below the scale /x = mw, the effective Hamiltonian for b — > sl + (r takes the form || 



n 



AG 



w 



-£-V tb V t :J2c^)0^), (4) 

* 1=1 



where we have used unitarity of the CKM matrix to remove V c bV* s dependence and have 
neglected the tiny V u bV* s terms. The operators Oi(ji) are 

Ox = {sL<xlnbLI3){cLf3l >1 C La ), 2 = {s L *l^ Lct ) {c Lp ^C Lp ) , (5) 

q=u,d,s,c,b q=u,d,s,c,b 
q=u,d,s,c,b q=u,d,s,c,b 

°i = T^^-^^Pr + rh s P L )b, 8 = -^s a T^ v (m b P R + m s P L )b p G a ^ , 

(? (? 
°v = TT-zS La 'fb L J'y ll £, O 10 = -r^s La Yb La h^, 

lorr lo7r 

where Pr s l = (li 7s)/2. In the following, we shall neglect the mass of the strange quark 
in 07.8- For our analysis, m s is not needed as a regulator for IR divergences, which are 
explicitly cut off by non-perturbative scales ~ Aqcd- In the shape function region, the m s 
dependence is small and was computed in Ref . j|6| . Non-perturbative sensitivity to m s shows 



up only at subleading power, while computable (m;?/m b Aq C d) jet-function corrections are 
numerically smaller than the Aqcd/^6 power corrections. 

At NLL order, one requires the NLL Wilson coefficient of Og and the LL coefficients of 
the other operators. For ©7,9,10 these are given by [JT0|, 



_16 X / _14 _16\ ^ — - 

C™> R (») = r »C 7 {M w ) + -[r » - r 23 J C 8 (M W ) + J> r^, (6) 

i=i 

Cr R M = P Nm V) + FK /f l) - + PeWEM/M*,), 



CM = C 10 {M W ) 



sin 9 W 



sin 2 



where C>?(mw), C%(mw) and the Inami-Lim functions Y, Z, and E are obtained from 
matching at /i — mw, and are given in Appendix |A|. The /i-dependent factors include JT^, [TT 



Pr R M = ( -0-1875 + J> r^' 1 ) + 1-2468 + (pf DR + s t r^) 



i=l I i=l 



Pb(ji) = 0.1405 + ft ^o ai_1 , r = -^L . (7) 

a s {m w ) 

The numbers t iy a i; pf DR , s», ^ that appear here are listed in Appendix [A]. Results for 
the running coefficients of the four-quark operators, Cx-&{n), can be found in Ref. JlIJ. We 
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have modified the standard notation slightly (e.g. 7"o(a0) t° conform with additional stages 
of the RG evolution discussed in sections [11 B| and [11 E[ Contributions beyond NLL will be 
mentioned below. 

At a scale \x ~ m;,, we need to match 6 — > s£ + £~ matrix elements of Tiw on to matrix 
elements of operators in SCET with a power expansion in the small parameter A, where A 2 = 
Aqcd / m b- For convenience, we refer to the resulting four-fermion scalar operators in SCET 
as "currents" and use the notation J«. In SCET we also need the effective Lagrangians. The 
heavy quark in the initial state is matched on to an HQET field h v , and the light energetic 
strange quark is matched on to a collinear field £ n . For the leading-order analysis in A/m^ 
we need only the lowest-order terms, 



n 



GpOi 



V2< 



71 



<y»K) j< 



(0) 



~ ^HQET ' ''-SCET > 



,(0) 



(8) 



where jf^ is the LO operator and the quark contributions to the HQET and SCET actions 
are 



C 



(0) 

HQET 



h v iv • D us h v 



(o) 

SCET 



1 



in-D c + i]£— r i]pk 
in- D r 



2 



(9) 



The covariant derivatives D us and D c involve ultrasoft and collinear gluons respectively, and 
we have made a field redefinition on the collinear fields to decouple the ultrasoft gluons at 
LO |[29||. For convenience, we define the objects 



T~C V Y h v , ipus Y q us , T) v 
Xn = W^ n , V C = W^D C W, igW c 

which contain ultrasoft and collinear Wilson lines, 



Y ] D US Y 

iw^[in-D c ,iD^W 



Y(x) = Pexpiig / ds n- A us (x+ns) ) 

J — oo 



(10) 



(11) 



and 



W(x) 



Pexpiig / ds n- A n (x+sn) ) 

J — oo 



(12) 



as well as the label operator V [£S 



To simplify the analysis we treat both m c and m& as hard scales and integrate out both 
charm and bottom loops at /x ~ ttv At leading order in SCET, the currents that we match 
on to are 



i=a,b,c 



i=a,b,c 



C 7j (s) 2m B (xn, p rf )fl H v ) (£ llM £) 



(13) 



J=a r 
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where the sum is over Dirac structures to be discussed below. The simple structure of these 
LO SCET operators is quite important to our analysis: for example, by power counting there 
are no four-quark operators that need to be included in SCET at this order. In Eq. flI5|) 
two auxiliary four-vectors appear, v M and n M . The B momentum, total momentum of the 
leptons, and jet momentum (sum of the four-momenta of all the hadrons in X s ) are 

p B = m B v^ , = p^+ + f t - , Px = n 'Px y + n-px — , (14) 

respectively. Here v 2 = 1 and and are light-like vectors, which satisfy n 2 = n 2 = and 
n-n — 2. The components of a vector can then be written as (p + ,p~ ,p±) = (n ■ p,n -p,p±). 
We use a frame in which = v ^ = and t> M = (n^ + n M )/2. Since px = thbV — q we have 

p\ = m 2 x = n-px n-px = m B + q 2 — m B {n-q+n-q) , q 2 = n-q n-q , 

n-px = rriB — n-q , n-px = ttib — n-q . (15) 

For later convenience we define the hadronic dimensionless variables 

2^£- _ n-px n-px q 2 , . 

xh = , Vh = > u h = , Vh = — ■ (16) 

m,B rriB ms m B 



In SCET the total partonic n ■ p momentum of the jet is a hard momentum ~ rrib and 
also appears in the SCET Wilson coefficients. At LO, n-p = (ml — q 2 )/mb and demanding 
that n-p is large means only that q 2 cannot be too close to m 2 . For example, neither q 2 w 
nor q 2 m m\j1 modifies the power counting for n-p. Thus, there is no requirement to impose 
a scaling that q 2 be small. For convenience, in the hard coefficients we write 



q 2 

C(n-p,m b , no, fi b ) -> C(s,m b ,p ,fx b ) , s = — , (17) 

m b 

since the partonic variable s is a more natural choice in b — > s£ + £~ and is equivalent at LO. 
For purposes of power counting in this paper we count s ~ A . We shall see in section [III Ej 
that varying s causes a very mild change in the coefficients. In Appendix [B] we briefly 
explore a different scenario, in which s ~ A 2 . A distinction between two matching scales hq 
and Hb is made in C in order to separate the decay rate into two /i-independent pieces, as 
displayed in Eq. (|3|). For power counting purposes, [Aq ~ \i b ~ m b and formally /^o > A*6- For 
numerical work one can take Hq = fi b . 

In Eq. ( pT3| ) we begin with a complete set of Dirac structures for the vector and tensor 
currents in SCET, namely 

ril = P R { r ,^,^}, r^P R %{i^,^,^, n ^}. (is) 

I n-v J q A v n-v n-v J 

These come with Wilson coefficients CQ a ,b,c and Ci a ,b,c,d respectively. This basis is over- 
complete for B — > A s £ + £~, but considering a redundant basis makes it easy to incorporate 
pre-existing perturbative calculations for the currents into our computations. Only the 
coefficients CV ai g a appear at tree level, but for heavy-to-light currents it is known that the 
other structures become relevant once perturbative corrections are included. For simplicity 
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of notation, we treat the 1/q 2 photon propagator in as part of the effective-theory 
operator. 4 

To reduce the basis in Eq. ( ^8|) further, we can use i) current conservation, q^i^^i = 0, 
ii) q^t^n^l = for massless leptons, iii) a reduction of the tensor Dirac structures into 
vector structures, since they are all contracted with q T . Constraint ii) allows us to eliminate 
Cioc- Taken together, constraints i) and iii) allow us to reduce the seven terms Cgi and C 7 i 
to two independent coefficients. For our new basis of operators we take 

4 0) = c 9 (xn, P P R rn v )(i^£) -c 7 ^^(xn, P PRi^ T n v )(I lfl £) 



+ ClOa {Xn^PRl^v) (tjvlst) + C 10b (Xn, p P R V^H v ) (t-f^d) , (19) 



and find that 



n n C 9b m B 2m B {C lc - C 7d ) + n-qC 9c 

L-9 — C 9a H Wb H Z > {^) 

2 n-q n-q— n-q 



C 7b n-q If -q 2 



2m 



B 



C 9c - n-qC 7c + n-qC 7d 



W — Wa « "j <^96 i Z 

2 Am B n-q— n-q 

ClOa = ClOa , 

r -r i 2ra " 9 r 

L>10b — ^10b H Z ^10c • 

Our Dirac structures for the Cg and C7 terms in Eq. ( |19D were deliberately chosen, in order 
to make results for the decay rates appear as much as possible like those in the local OPE. 
The fact that the basis of SCET operators for B — > X s £ + £~ involves only bilinear hadronic 
currents at LO means that in the leading-order factorization theorem we find the exact 
same non-perturbative shape function as for B — > X s 7 and B X u iv. This is immediately 
evident from the operator-based proof of factorization in Ref. p9"|| , for example. While the 
coefficients Cgi, C 7 i, Cioi in Eq. fll3|) are functions only of s = (n-q) (n-q) /m 2 ., the reduction 
of the basis of operators brings in additional kinematic dependence on n-q and n-q for the 
Cj's (which is also the case in analyzing exclusive dilepton decays [0). At tree level we have 
09,10 contributing to C 9a and C Wa , and a contribution from 7 with the photon producing 
an £ + £~ pair, which give 

Cg = Cg DR (fio) , C 7 = — C^ DR (/io) , CiOa = C10 , C\ob = . (21) 

m B 

Beyond tree level there will be C 7 dependence in Cg, and Cg dependence in C 7 . Eq. ( |2"Tp 
indicates that with our choice of basis the same short-distance dependence dominates in 
SCET: Cg ~ Cg, etc. We explore this further in Section |111 E] . In Eq. (E^) there is no 



distinction as to whether this matching is done at fi = //o or /i = fj, - The effective-theory 
operator in Eq. ( |T9"D was defined with a factor of m B pulled out so that the /z-dependent 
factors fribC^^ are contained in the coefficients C 7 . 



4 If we instead demand that the momentum q 2 be collinear in the n direction, with s ~ A 2 , then the SCET 
operator with a photon field strength should be kept, and will then be contracted with an operator with 
collinear leptons within SCET. In this case there will also be additional four-quark operators needed in 
the basis in Eq. (|T^). 
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FIG. 2: Graphs from H\y for matching on to SCET. 




FIG. 3: Graphs in SCET for the matching computation. 



At one-loop order, the full-theory diagrams needed for the matching are shown in Fig. |2| 
(plus wave-function renormalization, which is not shown). At this order the four-quark 
operators Ci_6 contribute through Fig. ||a. The one-loop graphs in SCET with the operators 
in Eq. (|19|) are shown in Fig. |3] (plus wave-function renormalization, which is not shown). 
There are no graphs with four-quark operators within SCET since we treat q 2 ~ A , so 
Fig. ^|a matches directly on to Cg. 

As discussed in the introduction, we perform a split- matching procedure from the full 
theory above m b on to SCET below mj, making use of two matching scales fi and /v 
Contributions from this stage of matching therefore take the form 

B(ji Q ,(i b ) = B 1 (ji Q )B 2 (ji b ). (22) 

Since Oig has no anomalous dimension above rn^ and there is a common universal anomalous 
dimension for all the operators in jf^ below m&, there is a well-defined prescription for 
carrying this out. We take all contributions that cause perturbative corrections to C Wa and 
Ciob to be at the scale so for this operator Bi(n ) = C±o, and at one-loop order i?2 (/•**>) 
includes a s (^b) In 2 (/•*&), Oi s (^b) m (/i&), and a s {^b) terms from matching the vertex diagram 
Fig. |||b and wave-function diagrams on to SCET. The analogous contributions from vertex 
diagrams for Cg and C-j are also matched at \x = \ib to determine their i?2(/i&)'s (for C-j the 
full-theory tensor current has a In /x that is matched at \x = /Jo) • The universality of the 
anomalous dimensions in SCET guarantees that this procedure remains well defined at any 
order in perturbation theory and can be organized into the product structure displayed in 
Eq. (^). For Cg and Cj there are additional non- vertex-like contributions that are matched 
on to -Bi(/io) at a scale /i > These include contributions from four-quark operators 
Oi_6 in the full theory, which will match on to Cg and C-j in SCET. 

The difference between the full-theory diagram in Fig. and the SCET graphs in 
Fig. |b, c is IR finite (where we must use the same IR regulator in both theories, as is 
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always the case for matching computations). In the UV the full-theory graph in Fig. |]b plus 
wave-function renormalization is /^-independent since the current is conserved. The graphs in 
SCET induce a \x dependence and an anomalous dimension for the effective-theory currents. 
These terms are matched at /i = We start with the basis in Eq. ([13]) and find 



ClOa(/^0) (l>b) 
Cl06,10c(/^Cb fJ'b) 



c 



10 



1 + 



w tt (S, fib) 



7T 



a s {fi b ) y 
Cio ojZas 



7T 



(23) 



with a constant /^-independent C w . The perturbative coefficients were computed in 
Ref. [27|, and setting n-p/rrib = (1 — s) we find 



u a 0, fib) 



21n 2 (l-s) + 2Li 2 (s) + ln(l-s) 



l-3s 



7l~ 



+ — + 6 
12 



2 In 2 ( !*) + 5 In ( - 4 ln(l - s) In ( B. 



«m 6 , 

-+ ? - ln(l - s 



(2a-l)(l-a) 



ln(l 



:i- s ) 



(24) 



For the matching on to Cg aj b, c in the basis in Eq. ([131) we have the same perturbative 
coefficients u a ^, c as for Cio a ,&, c , because only the leptonic current differs: 



C9a(/^o,/ i fe) — Cg nix (A t o) 

C%,9c(/^0; Mb) = GgVMo) 



1 + 



Ot s (Hb) 



71 



Ots{Hb) 



71 



b,c\ 



(25) 



However, for C 9i there are additional contributions, C™ lx (/i ), from the matching at /i = /a , 
which at one-loop order and C?(a°) includes Fig. |2|a: 



C 9 mix (^o) = C 9 NDR (/i ) + ~ (3C 3 + C A + 3C 5 + C 6 ) - s) (4C 3 + 4C 4 + 3C 5 + C 6 ) 



777 I 

+ /i(— , s) (3d + C 2 + 3C 3 + C 4 + 3C 5 + C 6 ) - -h(0, s) (C 3 + 3C 4 ) 
y m b ' 2 



+ 



a s (/f 



7T 



Oj r ,mix(l) 
°9 



(>o) 



(26) 



where all running coefficients on the RHS are C, = Cj(/io). We shall discuss the relation of 
C 9 mix to Cf in in the local OPE analysis ]TTJ, [TJ after Eq. ([33]). In Eq. (f§) the functions 
h(l, s), h(z, s), and h(0, s) for the b-quark, c-quark, and light-quark penguin loops are ||, [EL 



h(z, s) 



Mo 

9 v mft 



9 



2 
9 1 



0(1 -C) -m + ln 



1 + y^C 

1-VT^C 
4 



0(C- 1) 2arctan 



^^2^^ ln 0^ lnS+ 9- 



(27) 
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with ( = 4z 2 /s. Higher-order 0(a s ) corrections in Eq. (^) are denoted by the C™ 1X ^ term. 
An important class of these corrections from mixing can be determined from the NNLL 



analysis in Refs. [13], [14], [16]: 

^(^o) = C? DR ^8^9(s,/io) + Ci Ki^ 9 (s,fi ,rh c ) + C 2 H2^g{s, Ho, rh c ) . (28) 

To determine these terms one must be careful to separate out the factors in square brackets 
in Eq. (F25|). However we shall not attempt to include all NNLL terms consistently here. 
Contributions to C™ lx(1 ^ from the penguin coefficients C 3 _ 6 are unknown but expected to be 
small (at the ~ 1% level). 

Lastly, we turn to the results for C^. From the vertex graphs we have 



C7a{fJ>0, Hb) — C™ x (fio) 



a a (fj, b ) T 
1 + u a {s,(i b ) 

7T 

Ots{Hb) 



71 



(29) 



The ujJ perturbative corrections are again determined from the SCET matching in Ref. |p7 
which (switching to s) gives 



ul(3,fi, b ) = -^[21n 2 (l- S ) + 2Li 2 ( S ) + ln(l- S )( 



2-4s 



7T 

+ — + 6 

12 



+ 2 In 2 (J±) + 5 In ( - 4 ln(l - s) In ( ^ 

\m b / \m b / \m b 



a; 6 T ( S )= W J( S )=0. 



-2(1 - s) ln(l - s) 



(30) 



Additional contributions from other diagrams are matched at the scale /io into C™ lx (/xo)- 
Note that, unlike the vector currents, the tensor current for 7 gets renormalized for /i > m b 
and we must include the corresponding \B.(fj, /m b ) in C™ lx (/io), i-e. 



m b (fx ) 
m B 



i _ 2a a (/x ) ^ / /io \ 
3n \m b J . 



+ 



0^5(^0) £,mix(l)^ ^ 
7T 7 



where, much like in the case of C™ lx , we have 

C™ x(1 Vo) = Cg DR K 8 a (s, /i ) + Ci «a(s, /io, m c ) + C 2 k 2 (s, /i , m c 



(31) 



(32) 



and the results for k 8 ^ 7 (s, /i ), Ki-y^s, /io, ?^c), and ^2^7(5, /i , m c ) can be found in Ref. [^8 



Contributions to C 7 nix( ^ 1 ' ) from the penguin coefficients C%-q can be found in Ref. |49 



Using Eq. (|20|), n-qn-q/m 2 B = y H , and n-q/m B = 1 — we can use the above results 
to give the final coefficients for our basis of operators with the minimal number of Dirac 
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structures, namely 



c 9 = c 9 mix M(i + ^ 

7T 



[S, fl b ) + ~UJ b (S) + — 

2 (l-u H ) 2 -y H 



n (l-u H y-y H (l-u H )_ 



7T 



7T 



VhujJ{s)-{1-uh) 2 ujJ{s) 



(l-u H ) 2 -y H 



~ y H u^(s) y H (l-u H )uV(s) 
Ml-u H ) 2[(l-u H ) 2 -y H ] 



ClOb — Cio 



a s (Vb 



7T 



2(l-U H ) 2 y 

T\ ^2 U c (S) 

y l-u H ) 2 -y H 



(33) 



where the terms have the structure of a sum over products Bi(fi )B2(ftb), as desired. 

In using the results in Eq. (|33| ) one can choose to work to different orders in the /io- and 
/it-dependent terms, as shown in Eq. (|J). For the fig dependence, C™ lx (/io) and C™ lx (/io) 
include terms from matching at my/ and running to m&, as well as matching contributions 
at mb that cancel the fig dependence from the other pieces. Thus, these coefficients have 
only a small residual fig dependence, which is canceled at higher orders, just as in the 
local OPE. The Cj coefficients depend on fi b , both through a s (/x&) and through explicit /j 6 
dependence in u;J and lu^. The In /i 6 dependence in and cuj is identical, as expected from 
the known independence of the anomalous dimension on the Dirac structure in SCET. The 
lib dependence in Ci(ftb, /-to) is universal, and will cancel against the universal fib dependence 
in the jet and shape functions, which they multiply in the decay rates. We consider the 
phenomenological organization of the perturbative series for fig and fib terms in turn. 

First consider the /io terms. Because of mixing, the sizes of contributions to C<^ DR are 



comparable at LL and NLL orders |10], [11]], so a reasonable first approximation is to take 
the NLL result (just as for the local OPE decay rate). This entails dropping the 0(a s ) 
matching corrections C™ 1X ^ and C™^\ and running Cg at NLL order with Cj at LL order. 
As an improved approximation, we would then adopt the operationally well-defined NNLL 
approach [13[] of running both Cg and Cj to NLL order and keeping the 0(a s ) matching 
corrections at m^. 5 

Below mb there are Sudakov logarithms. For the fib dependence, the RG evolution in 
SCET sums these double-logarithmic series. As a first approximation we could take the LL 
and NLL running in UnifJ-b, Hi) and Us(fii, fi\) in Eq. (|j), while using tree-level matching for 
B 2 (fib) and J(fii). This is consistent because the NLL running is equivalent to LL running 
in a single-log resummation. As a second approximation we could then take NNLL running 
in both terms and include one-loop matching for both B 2 {fib) and J(fii). However since 
the scales ml 3> m b A ^> 1 GeV 2 are not as well separated as m\y ^> mf, we could instead 



5 We assume that matching at the high scale, mw, is always done at the order appropriate to the running 
of Uw(p>W,(M>) hi Eq. (||). 
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consider the second approximation to include the one- loop matching for i?2 (/•*&) and J(fii) 
with NLL running, but without including the full NNLL running (for which parts remain 
unknown) . 

Our procedure for split matching above was based on the non-renormalization of O\o in 
QCD. It can also be thought of as matching in two steps. First one matches at fio on to the 
scale-invariant operators 

J<°) = C 9 mix (sP R Yb)(h/) + C w (sP R rb)(I^ 5 £) 

to determine the coefficients C™g x . These coefficients are fio independent at the order in 
perturbation theory to which the matching is done. Secondly, the operators in Eq. ([53]) are 
matched on to the SCET currents in Eq. (|T9|) at the scale fib to determine the coefficients C7, 
C9, Cioa.fe- In Eq. (|4|) the operators for C™ lx and C\o are conserved, but the tensor current 
has an anomalous dimension, and so we take fi = as a reference point for matching on to 
a scale-invariant operator. This choice corresponds to the In mj factor in Eq. (|3l|) for C™ lx . 
A different choice will affect the division of a s (fi ) or a s (fib) terms. Note that Eq. (|34"D 
should be thought of only as an auxiliary step to facilitate the split matching; there is no 
sense in which the running of the tensor current is relevant by itself. In general the split- 
matching procedure could be carried out in a manner that gives different constant terms at 
a given order, but any such ambiguity will cancel order by order in Cj and C$ (and explicitly 
if fi = fi b ). 

Finally, note that our u a differs from the result for a; identified in Ref. for the 
partonic semileptonic decay rate when using the local OPE, 



OPE _ _1 

^semi g 



2 H s) K1-.) + 4Li 2 ( s ) + M l_.) (f±|) + K.) 



(5+9s-6s 2 ) 2tt 2 



2(l-s)(l + 2s) 



(35) 



Here w° m f contains both vertex and bremsstrahlung contributions evaluated in the full 
theory. Grouping these contributions with the Wilson coefficient for Oq gives 



&r\fi) = cr(fi) + p NDR (fi) ^uz^ ( 30 ) 



7T 



which is Cg S in the notation in Ref. ||10| . At LO, the restricted phase space in the shape 
function region causes bremsstrahlung to contribute only to the jet and shape functions, and 
not at the scale fi ~ m&. The shape function and jet function also modify the contributions 
from the vertex graphs. Thus, instead of cu^f the final results in the shape function region 
are given by our toj and ojJ factors appearing in C 9i and C 7i . Consequently, the main 
difference is in the terms we match at fi = fib, while the terms matched at fi — fio that 
appear in C™ lx and C™ lx are identical to terms appearing in the local OPE analysis. 



B. RG Evolution Between fib and /i; 



The running of the Wilson coefficients in SCET from the scale fif ~ m 2 to fif ~ m{,AqcD 



involves double Sudakov logarithms and was derived in Refs. [26, U% at NLL order. The 
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SCET running is independent of the Dirac structure of the currents, which is a reflection of 
the spin symmetry structure of the current. We briefly outline a short argument for why 
this is true to all orders in perturbation theory The leading-order currents in SCET have 
the structure 



J = {l n W) p Y{Y^h v ) 



(37) 



and we wish to see that their anomalous dimension is independent of V. The anomalous di- 
mensions are computed from the UV structure of SCET loop diagrams, with the Lagrangians 
in Eq. (Bp. Soft gluon loops involve contractions between the Wilson line and the h v and 
do not change the Dirac structure. Next consider the collinear loops. The attachment of a 
gluon from the Wilson line W to the collinear quark gives a factor of a projection matrix, 
which can be pushed through 7j_'s to give £ n f^/4 = £ n . Thus it does not modify the Dirac 
structure, so only insertions from the i$rl/(in- D c )ip^r term are of concern. These terms 
give structures of the form ■Un"' ) 7^_ 7^ 2 • • ■ r )^ k Tuh \ where all fii indices are contracted with 
each other. Using {7^,7^} = 2g^ and 7^7^ = d — 2, we can reduce this product to terms 
with zero 7j_'s since all vector indices are contracted. Hence all diagrams reduce to having 
the Dirac structure that was present at tree level, u^Tu^. 

Thus, all the LO coefficients obey the same homogeneous anomalous dimension equation, 



djj, 



r cusp (a s ) In (jt) + 7(a s ) d(^) 



(38) 



-r 



rn, P K<y s ) In ( — ) + \i{a a ) + r cusp (a s ) In ( — 



This must be integrated together with the beta function (3 = fj,d/dfia s (fj,) to solve for Uh in 

Ci(fii) = y/U H (fJ>i,fJ>b) Ci(fJ. b ) • (39) 

In the second line of Eq. ( [3~g|) we used the fact that V gives the total partonic n-p momentum 
of the jet X s in the B — > X s £ + £~ matrix element, and we introduced artificial dependence 
on the matching scale fib in order to make the n-p dependence appear in a small logarithm. 
Here n-p — — n-q. We write 



-\cusp 



n=0 



cusp / ^ s 

An 



n+l 



7 



E 

n=0 



In 



47T 



n+l 



P 



00 



n=0 



4:71 



n+l 



(40) 



At NLL order we need f% = HC A /3 - 2n//3, l3 x = 34C 2 A /3 - 10C A ^//3 - 2C F n f and 



pcusp = 



ri 



cusp 



8C pB = 8Cp 



67 tt 2 x 5 

Ca\ — n f 

1 18 6 / 9 i 



70 = -5Ci 



(41) 



50 



where Ca = 3 and Cf = 4/3 for SU(3). For the number of active flavors we take nf — 4 since 
we're running below m&. The cusp anomalous dimension r^ usp was computed in Ref. 
and the result for Y™ sp was recently found in Ref. [51 
order has been considered in Refs. pHI, I 



RG evolution in SCET at NNLL 



I For the NNLL result one needs r™ sp , ji, and j3 2 - 
For 7x an independent calculation does not exist, but a conjecture for its value was given 
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in Ref. [[44] based on the structure of the three- loop splitting function [51|. For the sake of 
clarity we stick to NLL order here. The result is 



U H (Hi,Hb) = exp 
where the independent variable is Hi and 



+ 2g 1 (r 1 ,n-p) 



2tt 



(42) 



a s (n b ) 2tt + f3 a s (n b )\n(ni/jd b ) 



with 



9o(n) 
gx{r u fi-p) 



— 1 + In T\ 



(43) 



(44) 



1 — r\ + ri lnri — - In r\ 



+ ■ 



Pa 



2 V fi b 



2C F B 



ri — 1 — In r\ 



This is the form for the universal running of the LO SCET currents found in Ref. p7 
Switching to a s (fii) as the independent variable, with r\ = a s (/ii)/a s (Hb), gives 



lllri 



/'6 



exp 



2ffo(r 1 ) 



(ri) 



(45) 



where go( r i) is as in Eq. ( [44]) and 



5C C 

,/, ( / , ) - — In 2 n + -/ In n + -J (25/5 - ft) (l - n + In n) 



(46) 



This form of the evolution with a s (n) as the variable was used in Ref. p4] , and is also the 
one we adopt here. The decay rate is computed from a time-ordered product of currents 
and so at the intermediate scale /if ~ m b A will involve products 



Ci(/ij, Ho) C s {ni, Ho) = U H (fii, Hb) Ci(/2 b , ho) Cji/Jb, Ho) 
explaining why we used a notation with \f\Ju in Eq. (^) . 



(47) 



C. Hadronic Tensor and Decay Rates 

In the last two sections we constructed the required basis of SCET current operators with 
matching at Ho ~ ~ m b anc ^ evolution to Hi ~ m b A. At the scale Hi we take time-ordered 
products of the SCET currents and compute the decay rates using the optical theorem. In 
this section we discuss the tensor decomposition of the time-ordered products and results 
for differential decay rates. 
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In order to simplify the computation of decay rates it is useful to write the sum of hadronic 
operators as a sum of left-handed and right-handed terms since for massless leptons we have 
only LL or RR contributions |30|. Doing this for our current, we have 

J$ = [Cg ~ Cioj {Xn %P L U v ) (£-fP L t) + [C 9 + C Wa ] (Xn l^L H v ) (I^Pr t) 



+c wb (xn v,p r n v ) (f y 75 1) - c 7 ^f- 

(Jl^L^ + Jr^L r ) , 



where 



J L{R) 



I = e-fp L t, L R = £YP R i, 

2m B ^ j 



(48) 
(49) 



XuPr 



7 



TCwbV 11 



Thus, the inclusive decay rate for B — > X s £ + £ is proportional to (W^ V L L V + W^ U L R U ), 
where the leptonic parts L^ R s and hadronic parts Wj% R s are given by 



[ 1 l(R)(p+)i'*Il(r)(P-)] [Il(r)(P-)YIl(r)(p+)] 



spin 



2 [p% p v _ +ptp u + - <T P+-P-T ie^ P+a p_ p ] 



(50) 



and 

W HR) 



^^(2nmp B - q -p x )(B\J^^X)(X\J^\B) 
B x 



Here, we use relativistic normalization for the \B) states. For convenience, we define pro- 
jection tensors Pf v so that 



MR) 



MR) 



(51) 



MR) 



yyL(R) _ piiu yyL{R) _ 



Hi/ 



(52) 



They are 



n 



fJIV 



r 2 



r A 



1 p Q 2 v ^ vV + Q^ff ~ v-q{v ,1 q v + v v q> 1 ) 
~2 9 + 2[q 2 - (vq) 2 } 



3q 2 P£_ + gV" - q"q v 
[q 2 - (v-q) 2 } 

g» v - v»v v + 
[q 2 - (v-q) 2 } 



pfiv 



-ie^qaVp 
2[q 2 -(vq) 2 } 



r 5 



2v-q 



(53) 



The optical theorem relates W^u to the forward-scattering amplitude defined as 



2m 



B 



d A xe- iq x (B\TJ^(x)J^(0)\B) 



-g^uTi + VpV v T% + ie^ a f]V a q p Tz + q^q v T^ + (v^q u + v v q^)T^ , 



(54) 



18 



with an analogous definition for Tf^, giving 



W 



ImT 1 



W. 



R 



71 



Im T; 

7T 



R 



(55) 



Contracting the lepton tensor L^ R ^ with and neglecting the mass of the leptons 

give the differential decay rate 



d 3 r 



96 



dq 2 dE_dE^ 



m 



B 



q 2 W x + (2E_E + - q 2 /2)W 2 + q 2 (E_ - E + )W 3 6{AE_E + - q 2 ) 



(56) 



where E± = v-p±,W\ = Wf+Wf 1 , W 2 = W 2 L + W 2 R , W 3 = W 3 L -W^ and the normalization 
factor is 



G 2 F m B a 1 



192vr 3 16tt 2 



\V tb V ts \ 



* |2 



(57) 



The Wi are functions of q 2 and v ■ q = v ■ (p + Another quantity of interest is the 

forward-backward asymmetry in the variable 



cos# 



v ■ p_ — v ■ p + 



where 9 is the angle between the B and i + in the CM frame of the 



(58) 



pair: 



d 2 A 



FB 



dv ■ q dq 1 



= / d(cos6) 



sign(cosi 



d 3 T 



48 q 2 



_ i T dv-q dq 2 d cos 9 

In terms of the dimensionless variables 

2E e - _ n-px 



ni 



[(v-q) 2 -q 2 ]W 3 . 



B 



X H 



m B 



Vh 



m B 



u H 



n-px 
m B 



(59) 



(60) 



the triply differential decay rate is 
1 d 3 T 



T dxn dy H du 



- = 2Am B (y H -u H ){(l-u H )(l-y H )W l + ^( y l-x H -u H )(x H +y H -l)W 2 



+ r ^(l-u H ){l-y H )(2x H +u H +y H -2)W 3 } 



(61) 



where Wi = Wi(uH,y~H)- F° r a strict SCET expansion we want n-p x ^ n-p x i.e. % y H . 
However, it is useful to keep the full dependence on the phase-space prefactors rather than 
expanding them, because it is then simpler to make contact with the total rate in the local 
OPE, as emphasized recently in Refs. |53|, [54j] , and so we keep these factors here. We shall 
also keep the formally subleading kinematic prefactors in our hard functions rather than 
expanding them as we did in Ref. [j38~| . Other variables of interest include the dilepton and 
hadronic invariant masses, 



Vh 



Q 



rn 



m 



2 ' 



Sh 



in 



x_ 

2 ' 



(62) 
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where 



sh = u H y H 



'l-u H ){l-y H ) 



(63) 



so that [y H > uh] 



1 - Vh + s H ± \f{l-y H + s H ) 2 - 4s H 



A few interesting doubly differential spectra are 
1 d 2 T 



T dy H du H 

i d 2 r 

T dy H ds H 
1 d 2 T 



24m B {y H -u H ) 2 {{l-u H )(l-y H )W 1 + l(y H - MH ) 2 ^ 2 } , 
2m B ^(l-y H + s H ) 2 -As H [l2y H W 1 + [(l-^ + ^) 2 -4^] W 2 } 



2rriB 



T dy H du H O—Uh)'' 



■[(l-u H ) 2 -y H ] \l2y H W 1 + 



[l-u H ) 2 -y H 



W 2 



d 2 T _ 2m B {s H -u 2 H ) 2 f 



To dsu dun 



it 



H 



0—uh) 

12u h {1-u h ){u h -sh)W 1 + (s H -u 2 H ) 2 W 2 



} 



For doubly differential forward-backward asymmetries we find 



d 2 A 



FB 
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u 
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(64) 



(65) 



(66) 



D. LO Matrix Elements in SCET 

At lowest order in the A/rrib expansion, the only time-ordered product consists of two 
lowest-order currents as shown in Fig. [|. The factorization of hard contributions into 
the SCET Wilson coefficients and the decoupling of soft and collinear gluons at lowest order 
are identical to the steps for B — > X s 7 and B — > X u £u, and directly give the factorization 



theorem for these time-ordered products [29]. The SCET result agrees with the factorization 



theorem of Korchemsky and Sterman |3"Tf . However, the structure of a s (\/mbA) and a s (rrib) 
corrections differs from the parton-model rate, as mentioned in Refs. [[33], |34|]. Beyond 
lowest order in a s (rrib) the kinematic dependences also differ, as mentioned in Ref. |[38|| . For 
B — > X u £u, the final triply differential rate with perturbative corrections at 0(a s ) can be 



found in Refs. 33, Bl. 



The factorization and use of the optical theorem is carried out at the scale /j, — //j, and 
we expand W { = wf ] + W*> 2) + . . . in powers of A = (Agco/m^) 1 / 2 (with no linear term). 



20 




FIG. 4: Time-ordered product for the leading-order factorization theorem. 
For B — > XJ + l~ we have bilinear hadronic current operators in SCET in Eq. flT9"| ) and so, 



as is the case for B — > X u £u, we find 



(0) 



hi(p$,Px>K) / dk + J^(p-,k+, f i i )f^(k + +A-p+, fH ) 



(67) 



This result is important, since it states that the same shape function appears in B — > 
X s £ + £~ as appears in B — > A, 7 and B — > A u ^z/. This formula relies on the power counting 
s ~ y H ~ A° that we adopted (and would not be true for the counting s ~ A 2 discussed 
in Appendix [Bp. At tree level the structure of this factorization theorem is illustrated by 
Fig. [|. The hard coefficients here are 
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P.^ - ^Tr 
3 4 




^3 



(68) 



with P v = (1 + ^)/2 and T = 7°r t 7°. In Eq. (|7D we have the same leading-order shape 
function as in B — > A s 7 and P — > A n £z/, namely 



/ (0) (^) 



1 Aiar 



2 J 4vr 

- (B v \h v 5(£ + -in-D)h v \B v ) 



(5„|W„(x)W„(0)|5„) 



(69) 



where x^ = n-xn^/2. This function was first discussed in Ref. |[24|| . The jet function i 
defined by J^(p~, k+) = (-1/tt) Im (Jfe+) x%+ - fc+), where 



i(0|T[x^,(0)^-Xn l( ,/(ar)]|0> = 6(u-u')5 2 (x ± )6(x 
and is known at one-loop order [33, 34 1, namely 

4 



■dk+ 
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-ik+x~/2 j(V)(h+ 



(70) 



¥x{ 5{z) 

a s (ni)C F 



1 + a ^° F (21n> £P£. - 31n^ + 7 - tt 2 



47T 



47T 

4 In z 



+ (41n^-3 
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(71) 



Ml-*), 
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where z = k + Despite appearances, only the combination zp\ appears in apart 
from the 6(1 — z). This last ^-function is induced by the soft function and, when one takes the 
imaginary part of the full time-ordered product, affects the complex structure. Therefore, 
we include it in our definition of J^°\p~ , k + ). 



E. RG Evolution Between fi\ and /x. 



The function cannot be computed in perturbation theory and must therefore be 
extracted from data. This same function appears at LO in the B — > X s 7, B — > X u lv 
and B — > X s £ + £~ decay rates. In practice, a model for /(°) is written down with a few 
parameters, which are fitted to the data. The support of /^(A — r + ) is — oo to A since 
r + G [0, oo). It is often convenient to switch variables to f^°\r + ) = f(°\A — r + ) which has 
support from to oo, although we shall keep using /(°) here. A typical three-parameter 



model is p6 



/ W (A-r+ Ai A ) = /<°>(r+ Wi) 



a b (r+) 



+ \b-l 



r(6) L b 



exp 



-ar 



L 



6(r 



(72) 



where a, b are dimensionless and L ~ Aqcd- These parameters can be fitted to the B — > X s 7 
photon spectrum and the function /W can then be used elsewhere. The most natural scale 
to fix this model at is \i — ~ 1 GeV, at which it contains no large logarithms. The result 



of evolving the shape function to the intermediate scale is then |34 



' V' / (0) (A-r+',/, A ) 
r(r?)7 fil(r+-r+>y-v 



(73) 



(The structure of this result also applies at higher orders in RG-improved perturbation 
theory [fH}| , and at one-loop order a similar structure was considered earlier, in Ref. [|55|j.) 
At NLL order 
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(74) 



Here, r 2 = a s (/iA)/a s (/ij), TQ Uhp and r^ usp are the same as in Section [11 B| and 70 = —2Cf- 
For numerical integration this can be rewritten in the form 



/(°)(A-r+^) = e Vsi ^ A) 



r ( i + 1 ,)v£) , / 1<i(/, ° ,(X - r+(1 - il/ " ) '" A) - (75) 
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III. B XJ+£~ SPECTRA IN THE SHAPE FUNCTION REGION 



A. Triply Differential Spectrum 

At lowest order in the power expansion, Eqs. ( |55| ) and fl5"7|) give the result 



hipx^m^Hi) / dk + J ( - \p-,k + , tM )f ( - \k + +A-p x , f M i ), (76) 



(77) 



where RG evolution from the hard scale to the intermediate scale gives 

and the results at \i = \i b are determined from the traces in Eq. (|68|): 
1 / 1 1 2 I |2 \ 2Re\C 7 C*} 2\C 7 \" 

h( Px , Px ,» b ) = -{\c 9 \ +\c 10a \)+ +— - 

2(l-«if)/|„ |2 , -. ,> . r ,\ C 



(78) 



h 2(PxiPxiVb) 



-4Re[C 10a C 7 * 



(|C 9 r+|C 10a | 2 +Re[C 10a C i g) + 
2 Re[Ci 0a C 9 *] 
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8|C 7 | 



2 {i-yH){yH-u H ) 



m B (l-y H )(y H -u H ) m B (y H -u H )' 



Here Cj = Ci(p x ,p x , fi b , /io, ^&), so these hard coefficients also depend on m?, and have 
residual /io scale dependence. Explicit formulae are given in Eq. fl33l) . For convenience we 
define 



Px 



F(°\p x ,p x ) = U H ( f i t ,fi b ) / dk + J^\p-,k + ^ t ) f^(k++A-p x ,^) 

Jo 

= p$U H (m,Hb) [ dzJ^\p- ) zp x ^ l )f^rR-p x (l-z)^ l ). (79) 
Jo 

where p x = p~ + A. In terms of this function, 

W} 0) =h l (p x ,p x , f i b ) F(°\p x ,p x ). 

We find that to NLL order 



(80) 



F {0) (p x , Px ) = U H ( f i i ,ii b )f^(A-p x ,fi i ) 



(81) 
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Note that, until we include the a s corrections from the jet function, is independent of 
p x , so that all of this dependence is in the hi(p x ,p x , fi b ) functions. 
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Now, the triply differential decay rate in Eq. (|6"1~D becomes 

tt~ ~~j 1— — i — = 24m B (y H -u H ){ (l-u H )(l-y H )hi + -{l-x H -u H ){x H +y H -l)h 2 

1 o d%H dy H dun •> 2 

+ ~Y ^-UH)(l-yH)( 2x H + UH+yH- 2 )h^F {0 \m B u H ,m B y H ) , 

(82) 

with ^1,2,3 from Eq. (f7g|) . As a check on this result, one can make the substitutions 

C-9a = —ClOa = 1/2 , C 7 = Ciob = , (83) 

after which the foi and /i 2 terms in Eq. ([82]) agree with terms in the leading-order shape- 
function spectrum for B — > X u £z/ fl33| , |56| . The /i3 term for B — > X s ££ was the difference of 
products of left- and right-handed currents and so should not agree in this limit. 



B. d 2 T/dq 2 dm x Spectrum with q 2 and mx Cuts 

Next we discuss doubly differential rates and forward-backward asymmetries. For 
d 2 T /dq 2 dm 2 x the rate is obtained from Eq. ( |8"2"D by integrating over x B and changing vari- 
ables. In terms of dimensionless variables y B = <f /Trig and sh = rn 2 x /m 2 B we have 



1 d 2 Y 



H ys (y H , s H ) m B F (0) (m B u H (y H , s H ), m B y H (y H , s H fj , (84) 



T dy H ds H 

K ys {y H , s H ) m B F {0) (m B u H {yH, s H ),m B y H (y H , s H , 



1 d 2 A FB 
T dy H ds H 

where 

H ys (y H , s H ) = 2y/(l-y H +s H ) 2 -As H [uynh + [(1 -y H + s H f ~^s H ] h 2 ) , (85) 
K ys (y H ,s H ) = QyH[(i-yH+s H ) 2 -4s H ]h 3 

and we need to substitute ^1,2,3 from Eq. ([T8|) and u B (yH, sh) and y H (yH, s B ), as given in 
Eq. (p^p. When one takes experimental cuts on q 2 and m 2 x . 



yf <y H < yT x , < s H < s° H , (86) 
the limits on the doubly differential rate and forward-backward asymmetry in Eq. (^4j) are 



1) yf n <VH< yT X , 0<s H < mm{s° H , (l-^) 2 } , 

2) 0<s H <s° H , yT<y H <min{y^\{l-^) 2 }, (87) 

depending on the desired order of integration. 
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C. d 2 T/dm 2 x dp x Spectrum with q 2 and mx Cuts 

The hadronic invariant-mass spectrum and forward-backward asymmetry can be obtained 
by integrating the doubly differential spectra 

1 d 2 T 



H s (s H , u H ) m B F i0) (m B u H , m B —) 



To dsndu 
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D. d 2 T/dq 2 dp x Spectrum with q 2 and mx Cuts 

From Eqs. (|65|) and the above results, we can obtain the dilepton invariant-mass spectrum 
and forward-backward asymmetry, for example by integrating the doubly differential spectra 

^ = H«(y Hj u H ) m B F<f» (m B u H , m B izMZ^) (91) 

1 d 2 A FB (0) / l-y H -u H \ 
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1 o dyffduu V 1 — u B ' 

over u B - Here 

H y (yH,u H ) = m ~^ H)2 ~l H? \y H [(l-u H ) 2 +2yH] (|C 9 | 2 + |C 10a | 2 ) (92) 

+ [8(l-u H ) 2 +^y H ] \C 7 \ 2 + 12y H (l-u H )Re[C 7 C 9 *] 
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and the limits of integration with cuts are 



VT <VH< VT X , < u H < min jl- ^ , 

The opposite order of integration is also useful: 

< u H < 1 , yi(u H ) < Vh < y-i{uH) , 

(l-u H )(u H -s° H 



(93) 



yi{u H ) = max y% 



uh 



y 2 {u H ) =min|^ ax , (l-u H ) 



(94) 



mix /^mix 
9 j °7 ) 



The doubly differential rate can also be expressed in terms of the coefficients C, 
and Ciq. This is one step closer to the short- distance coefficients Cq, Cj, and C\q of Hw, 
which we wish to measure in order to test the Standard Model predictions for the corre- 
sponding FCNC interactions. Substituting Eq. ([33|) into Eq. (|92|) gives 



H y (y H ,u H ) 



i[(l-u H ) 2 -y H ]' 
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(95) 



where s = q 2 jm\ and 
n A = 14 



us 

IX 



(96) 



a s (fi b 
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7T 



u b {s,n b ) 



IX 

ot s (fi b ) 



7X 



2(l-u H y 
wZ(s, fi b ) - ujJ{s, fl b ) - UJ%{s, fl b ) 

UJ a (s, fi b ) - UJ C (s, fl b ) — r^— U b (S, fl b ) 



(2Q A n D + n B n c ) /3 . 

This is the form that turned out to be the most useful for the analysis in Ref. [30 



E. Numerical Analysis of Wilson Coefficients 

As shown in Fig. [3], for the small-g 2 window (q 2 < 6 GeV 2 ) we have p\ <C p~x- Generically, 
the hard contributions in C 9 , C 7 , and Cioa.iofc from our split-matching procedure depend on 
the variable q 2 . In Fig. |5] we plot the q 2 dependence of the real part of the coefficients 
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FIG. 5: Comparison of the real part of Wilson coefficients at //q = //& = 4.8 GeV with m c /m b = 
0.292, m&(//o) = 4.17GeV, and m b = 4.8 GeV. For Cg, C 7 , and Ciofe we take = 0. 

and see that there is in fact very little numerical change over the low-g 2 window. Here 
Re[C 9 ocal ] varies by ±1.5%, Re[C 9 mix ] by ±1%, and the real parts of {Cg, C 7 , C Wa , C 10b } by 
{±1%, ±5%, ±2%, ±3%}. The imaginary parts are either very small or also change by only 
a few percent over the low-g 2 window. The analytic formulae for the q 2 dependence mean 
that there is no problem keeping the exact dependence, but this does make it necessary 
to perform integrals over regions in q 2 numerically. A reasonable first approximation can 
actually be obtained by fixing a constant q 2 in the hard coefficients, while keeping the full 
q 2 dependence elsewhere. 

Since the coefficients change very little with q 2 we continue our numerical analysis by 
fixing q 2 = 3 GeV 2 . If we then take /io = Hb = n^b = 4.8 GeV, m^(/j ) = 4.17 GeV, 
m c /m b = 0.292 and p\ = we find that Eq. (|33"D gives 



C 9 = 0.826 C 9 mix + 0.097 C? 
C 7 = 0.823 C 7 mix + 0.001 C 9 n 



nax _ 9 



£<NDR 



/~<mix 

0.030 7 



£<NDR ' 



z^mix /^mix 

-0.239 ~^vtvd + 0.005 9 



Cj DR 



Cg DR 
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These numbers indicate that, despite the entanglement of C™^ in C 7i g due to a s (mb) cor- 
rections, numerically Cg is dominated by Cg and C-j is dominated by C7 in the Standard 
Model. 

For the coefficients at q 2 = 3.0 GeV 2 , with the other parameters as above, we have 



Cmix 
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Cg{u H = 0) 
C 9 (U H = 0.2) 
C\Qa 



4.487 + 0.046i, 
3.683 + 0.038i, 
3.663 + 0.038i 
-3.809, 



q mix 

C 7 (u H = 0) 
C 7 (u H = 0.2) 

C wb {u H = 0) 
C wb (u H = 0.2) 



-0.248, 

-0.198 + 6 x 10" 
-0.193 + 10 _4 z, 
0.214, 
0.237. 



(98) 
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The relevant range of p x in Fig. |l] gives < uh < 0.2. From the above numbers it is easy 
to see that the uh dependence of Cg, C7, and C106 is very mild over the range of interest. 
The perturbative a s corrections due to u( ,T reduce both Cg and C7 by 17% relative to C™ lx 
and C™ 1 * respectively, and C Wa by 15%. This can be seen both in Fig. | and in Eq. (p8|). 
when one notes that C\q = —4.480. Comparing with coefficients in the local OPE, we 
note that the co>° m f factor, which accounts for the difference between C 9 ocal and C^ 1X , is 
significantly smaller than the combination of a s corrections in the uj terms that shifts Cg 
from its lowest-order value. 

In quoting the above numbers, we have not varied the scales /io and yn>. The main point 
was to compare the size of the hard corrections in the shape function and local OPE regions, 
and to see how much deviation from C™g x they cause. The dependence on /xq for the C; is 



similar to that in the local OPE analysis at NLL [10, 11] and will be reduced by a similar 



amount when the full NNLL expressions are included in C™g x . The /i& dependence of the Cj 
is fairly strong because of the appearance of double logarithms, but it is canceled by the ^ 
dependence in the function F^°\ which contains the NLL jet and shape functions. 



IV. CONCLUSION 

In this paper we have performed a model-independent analysis of B — > X s £ + £~ decays 
with cuts giving the small-g 2 window and an mx cut to remove b — > c backgrounds. These 
cuts put us in the shape function region. We analyzed the rate for the formal counting 
with q 2 ~ A and m 2 x ~ A 2 and showed that the same universal shape function as in 
B — > X u tv and B — > X s 7 is the only non-perturbative input needed for these decays. We 
also developed a new effective-theory technique of split matching. Split matching between 
two effective theories is done not at a single scale /1, but rather at two nearby scales. For 
B — > X s £ + l~ this allowed us to decouple the perturbation-theory analysis above and below 
m&, which simplifies the organization of the a s contributions. 

In Section [TTT] we presented the leading-power triply differential spectrum and doubly dif- 
ferential forward-backward asymmetry with renormalization-group evolution and matching 
to 0(a s ). Above the scale mj, we restricted our analysis to include the standard NLL terms 
from the local OPE, but illustrated how terms from NNLL can be incorporated. Below m& 
we considered running to NLL and matching at one-loop (NNLL evolution will be straight- 
forward to incorporate if desired). We then computed several phenomenologically relevant 
doubly differential spectra with phase-space cuts on q 2 and m x (from which the singly dif- 
ferential spectra can be obtained by numerical integration) . In section |111[E we discussed the 



numerical size of our perturbative hard coefficients and compared them to the local OPE 
results. 

Our results for the doubly differential rate in Eqs. (|9l| ) and (p2[) , together with from 
Eq. (jSll), determine the shape-function-dependent rate for B — > X s l + l~ . Using as input 
a result for the non-perturbative shape function /(°) from a fit to the B — > X s ^ spectrum 
or from B — > X u iv gives a model-independent result for B — > X s i + i~ with phase-space 
cuts. A full investigation of the m^-cut dependence and phenomenology is carried out in a 
companion publication . An intriguing universality of the cut dependence is found, which 
makes the experimental extraction of short-distance Wilson coefficients in the presence of 
cuts much simpler. An extension of the analysis of this paper to include subleading shape- 



function effects will be presented in the near future [[57] . 
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APPENDIX A: WILSON COEFFICIENTS 



The coefficients and functions that appear in Eq. (H) are defined as follows fT0| 
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APPENDIX B: THE CASE OF COLLINEAR q 2 

In the body of the paper we used q 2 ~ A . We were free to choose this counting since 
the power counting for the leptonic variable q 2 does not affect the counting for in the 
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FIG. 6: Additional graphs in SCET for the matching computation for the case where q 2 ~ A 2 . 

shape function region. (The only restriction was not to have q 2 too close to m 2 .) However, 
we are free to consider other choices. In this appendix we consider how our analysis will 
change if we instead take q 2 ~ A 2 . With this scaling, new physical degrees of freedom are 
needed at leading order in SCET, making the analysis more complicated. In particular we 
must consider graphs with quark fields that are collinear to the collinear photon (or dilepton 
pair), since with this power counting we have (q ) 2 ^> q 2 . 

An example of a new nonzero graph is the one generated by four-quark operators within 
SCET, as shown in Fig. |], which involve these additional degrees of freedom. In this graph 
we have a light-quark loop of collinear-n fields that are collinear to the virtual photon. The 
presence of this type of diagram changes the hard matching at fib = tui,. It also means that 
we have a more complicated pattern of operator mixing within SCET, since divergences in 
the displayed diagram will cause an evolution for Cg, etc. Therefore, the running below m& 
will no longer be universal. In the presence of these diagrams the jet function will also no 
longer be given by a single bilinear operator, since it will also involve some contributions 
with a factorized matrix element of n-fields, which are also integrated out at p 2 ~ mbAqcB- 
Finally, the appearance of these additional degrees of freedom might also affect the number 
of non-perturbative shape functions that appear in the factorization theorem. It would be 
interesting to carry out a detailed analysis of this q 2 ~ A 2 case in the future. 

In B — > Xg'-f at lowest order, the analog of the graph in Fig. |6| vanishes at one-loop order, 



and this argument can be extended to include higher orders in a s |44]. This relies on the 
fact that here q 2 = and does not generate a scale. We find that the same reasoning does 
not apply for B — > X s ££ for parametrically small but finite q 2 . 

Finally, we comment on the possibility of penguin charm-loop effects. In our analysis we 
integrated out the charm loops at the same time as the bottom loops. This is reasonable 
when treating q 2 ~ A . One could also consider the case m 2 ~ m^A, which is also reasonable 
numerically. This type of power counting was considered for the simpler case of B — > X c iv 



decays with energetic X c in Ref. |58| and it would be interesting to extend this to B — > X s ££. 



We remark that the problematic region for B — > nn factorization theorems [59, BG, 61 



2[, which is near the charm threshold, q ~ 4m 2 , is not relevant for our analysis. The 



experimental cuts on q explicitly remove the known large contributions from this region. 
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